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摘 要 : 利用 ITS_LIVE 数 据 .Landsat 数 据 提 取 了 喀 喇 昆 
合 起 来 ,分 析 对 比 了 不 同 运动 类 型 冰川 表 
1989 一 2018 年 ， 研究 区 42 条 冰川 中 ,存在 稳定 冰川 16 条 前 进 冰 川 6 条 、 退 缩 冰 放 
云 动 一 般 原 理 ,可 在 
区 域 多 出 现在 跃 动 冰 川 主干 或 支部 ,流速 大 小 一 般 在 100.00 mea ' E, WEF 


动 与 表面 运动 特征 结 


条 。 空 间 上 ,流速 分 布 符合 冰川 运 


包谷 山 北 坡 和 2 条 冰川 的 表面 流速 。 将 末端 运 
流速 的 时 空 变 Lo 结果 表明 : (1) 
1 条 、 跃 动 冰川 19 
只 累 区 和 冰 舌 上 部 发 现 明显 的 高 值 区 域 ,这 些 

苏 盖 提 冰 川南 分 支 ， 


其 最 大 流速 可 达到 358.33 m- a, 2) 冰川 流速 分 布 与 地 形 要 素 的 关系 密切 。 流 速 在 海拔 4600~ 
5000 m 之 间 达 到 最 大 (54.55 ma'), 是 冰川 末端 流速 的 10 倍 以 上 ;分 布 在 坡度 0~5° 之 间 的 42 条 冰 


川 平均 流速 最 大 ,并 且 随 着 坡度 的 增 大 ,流速 逐 


流速 最 小 。(3) 稳定 冰川 流速 年 际 变化 较 稳 定 , 不 同年 份 相同 位 置 的 流速 值 较 
,年 际 流速 波动 均 较 大 ; 跃 动 冰川 不 同位 置 年 际 流速 变化 不 同 。 末 端 发 生前 进 


川 和 退缩 冰川 而 言 


逐渐 减 小 ;处 于 东 向 的 冰川 流速 最 大 ,处 于 西南 向 的 
一 致 ;而 对 于 前 进 冰 


的 跃 动 冰川 表面 流速 变化 不 大 ,但 末端 变化 速度 却 可 超过 100.00 ma^, E E TT 3& 8| 500.00 ma '。 


(4) KE Kr 
原因 。 

关 键 词 : 
文章 编号 : 


冰川 是 冰冻 圈 系 统 的 重要 组 成 部 分 ,同时 也 是 
气候 变化 指示 器 。 冰 川 变化 受气 候 变化 控制 ,并 
对 气候 变化 作出 快速 响应 , 主要 造成 冰川 积累 量 与 
消融 量 的 变化 而 引起 冰川 高 度 和 冰川 运动 的 变化 ， 
最 终 表 现在 冰川 面积 .长 度 等 发 生变 化 。 因 此 ,人 研 
究 冰川 变化 对 区 域 乃 至 全 球 性 气候 变化 具有 重要 
意义 。 冰 川 运动 是 冰川 的 一 个 重要 变化 特征 , 反 
映 了 冰川 在 内 外 环境 综合 作用 下 的 动态 变化 信 
息 “。 基 底 滑 动 与 底部 沉积 层 变 形 造成 的 运动 构 
成 了 冰川 的 主要 运动 形式 ”。 国 际 上 ,冰川 研究 起 
步 早 ,多 集中 在 极地 冰 盖 与 海洋 性 冰川 REK 
川 运动 观测 研究 起 步 较 晚 ,研究 主要 集中 在 天 山 ” 
祁连山 “中 喜马拉雅 山下 横断 山下 \ 喀 喇 昆 仑 
山王 ”等 地 区 部 分 冰川 ”“。 井 哲 帆 等 "利用 通 径 分 
析 方 法 发 现 , 冰 川 形态 .物质 平衡 .气候 条 件 等 对 冰 
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只 更 大 的 冰川 易 发 生 跃 动 。 气 候 、 雪 月 和 和 地形 也 是 影响 冰川 运动 速 


冰川 ; 喀 喇 昆仑 山北 坡 ; ITS_LIVE 数据 ; Landsat 数 据 ， 表 面 流 速 


度 变化 的 主要 


; 时 空 变化 


川 运动 速度 均 有 影响 , 强 弱 不 一 
度 对 运动 速度 影响 最 为 显著 。 
冰川 运动 特征 变化 与 冰川 物质 平衡 变化 存在 
一 定 的 因果 关系 ”。 在 全 球 气候 变 暧 背景 下 ,青藏 
高 原 冰 川 整体 呈 退 缩 趋 势 , 物 质 平衡 呈 负 。 但 由 于 
受 局 部 冰川 气候 条 件 影响 ,出 现 了 有 些 地 区 冰川 物 
质 平衡 呈正 的 现象 ,影响 了 冰川 运动 特征 ' ; 
Hewitt 等 ” 提出" 喀 喇 昆仑 异常 ”, 发 现 喀 喇 昆仑 
一 些 地 区 存在 冰川 停 讲 和 前 进 现象 。 duis H 
出 喀 喇 昆仑 山 冰 川 正 的 物质 
的 稳定 和 前 进 , 而 周边 区 域 出 现 冰川 后 退 和 减 薄 。 
许 艾 文 等 “利用 遥感 技术 对 喀 喇 昆仑 山 克 勒 青河 
流域 冰川 变化 进行 了 监测 与 研究 ,发 现 冰 川 退缩 速 
率 由 快 至 慢 ,并 存在 大 量 前 进 与 跃 动 冰川 ,造成 了 
冰川 面积 的 扩张 。 目 前 ,冰川 异常 运动 的 机 理 尚 不 


,但 冰川 面积 与 厚 
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李 海 娟 等 : 喀 喇 昆仑 山 乔 戈 里 峰 北 坡 主要 冰川 运动 特征 


清楚 ,但 影响 最 大 的 因素 是 受热 调节 而 不 是 水 文 控 
制 中 。Quincey 等 “指出 喀 喇 昆仑 山 有 跃 动 冰 川 存 
在 ,不 完全 符合 热力 学 和 水 文 控 制 , 最 大 的 差别 之 
处 在 于 冰川 基础 条 件 不 同 。Farinotti 等 指出 喀 喇 
昆仑 山 异常 现象 是 否 会 继续 持续 尚 不 确定 ,但 鉴于 
目前 对 未 来 气候 的 预测 所 预期 的 大 幅 升 温 , 这 种 异 
常 现 象 长 期 持续 下 去 似乎 不 大 可 能 。 

喀 喇 昆仑 山 冰 川 发 育 广泛 。 关 于 喀 喇 昆仑 山 
冰川 研究 也 较 多 ,其 北 坡 异常 的 扩张 现象 引起 了 学 
者 们 广泛 的 关注 。 不 少 学 者 研究 该 区 不 同 流域 的 
冰川 面积 .长 度 与 物质 平衡 状态 ,但 对 其 表面 运动 
研究 不 够 全 面 ” *”" , 且 未 进行 连续 观测 ,并 缺乏 对 
不 同 末端 变化 的 冰川 表面 运动 的 变化 特征 及 其 影 
响 因 素 的 监测 。 鉴 于 此 ,本 文通 过 人 研究 冰川 运动 与 
其 表面 地 形 之 间 的 关系 为 认识 该 地 冰川 异常 现象 
提供 资料 ,完善 该 区 冰川 运动 监测 。 


1 研究 区 概况 


喀 喇 昆仑 山北 坡 大 致 在 75°35’~77°00'E 和 
35°31’~36°?49'N( 图 1)。 据 中 国 第 二 次 冰川 编目 ， 
该 地 区 的 冰川 编码 为 “5Y654”, 流域 面积 为 8223 
km 。 喀 呵 昆 仑 山 主 山 迪 海拔 在 7000~8000 m 之 
间 中 。 该 区 现代 冰川 发 育 的 山脉 ,平均 高 度 均 在 
5730 m 以 上 ,主峰 乔 戈 里 峰 海 拔 为 8611 mm。 海拔 最 
高 的 冰川 为 乔 苞 里 冰川 ,高 程 为 8569 m, 其余 冰川 
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注 :图 中 数字 0-~41 代 表 冰 川 序号 。 
图 1 喀 喇 昆仑 山北 坡 冰 川 概况 图 


Fig. 1 Glacier profile of the north slope of Karakoram 


Mountains 


位 于 海拔 3979~8569 m ZZ [B] ,平衡 线 大 致 位 于 5$100-~ 
5400 m 之 间 ”。 与 周围 冰川 地 理性 质 比 较 ,该 区 冰 
川 由 复合 型 向 大 陆 型 过 渡 中 。 流 域内 音 苏 盖 提 冰 
川 是 中 国 面积 最 大 的 冰川 i。 

喀 喇 昆仑 山北 坡 气 候 受 亚洲 季风 影响 ,80% 的 
降水 是 由 夏季 降水 贡献 。 冬 季 受 西风 环流 控制 2/3 
的 降雪 发 生 在 高 海拔 处 下。 在 谷地 中 年 降水 量 不 
FE 100 mm, 但 随 海拔 升 高 ,降水 量 逐 渐 增 加 。 雪 线 
处 ,夏季 平均 气温 为 1.2~3.0% ,年 降水 量 在 500 mm 
以 上 '。 冰 川 物质 主要 来 源 于 雪崩, 雪崩 供给 了 冰 
J Bet, AS SCHR DE KI UE KU] .形态 选取 了 
42 条 冰川 作为 研究 区 (图 1)。 


2 数据 与 方法 


2.1 数据 来 源 

Landsat 数据 来 源 于 地 理 空 间 数 据 云 网 站 
(http:/Awww.gscloud.cn/) ,影像 行列 号 为 148/35 ,主要 
由 1989 一 2018 年 的 Landsat 5 TM Landsat 7 ETM+ 和 和 
Landsat 8 0LI 数 据 组 成 ( 表 1)。 研 究 所 采用 的 产品 
级 别 为 Level LIT 数 据 , 均 经 过 几何 校正 、 辐 射 校正 
和 基于 地 表 数 字模 型 的 地 形 校正 的 。 该 区 海拔 高 ， 
降雪 丰富 zs mA BUR KU us ,造成 解 译 困 
难 。 因 此 ,本 研究 尽 可 能 选取 夏 、 秋 季 云 筋 雪 餐 善 
较 少 的 影像 。 

ITS_LIVE 数据 (The inter-mission time series of land 
ice velocity and elevation ) 产 品 是 一 组 主要 冰川 的 年 
度 平 均 地 表 速 度 的 区 域 汇 编 ,时 间 跨 度 为 1985 一 
2018 年 。 下 载 地 址 https://its-live.jpl.nasa.gov。 它 是 
NASA“ 人 研究 环境 制作 地 球 系 统 数据 记录 2017" 项 目 
的 一 部 分 ,其 展示 的 地 表 速 度 是 使 用 Gardner 等 3 
描述 的 自动 裂 颖 特征 跟踪 处 理 链 的 方法 从 landsat 
4.5.7 和 8 图 像 获得 。 数 据 分 辨 率 为 240 m。 本 文 
选取 1990 一 2018 年 流速 数据 来 研究 冰川 表面 运动 。 

此 外 ,数字 高 程 模型 (DEM ) 数 据 是 来 源 于 地 理 
空间 数据 云 的 ASTER GCDEM 数 据 , 分 辨 率 30 m。 下 
载 地 址 (http://www.gscloud.cn/search)。 该 数据 主要 
用 来 提取 研究 区 地 形 参 数 ( 高 程 坡度 、 坡 癌 等 )。 
Sam 等 5 和 张 莎莎 等 已 研究 过 该 数据 的 可 靠 性 ， 
并 成 功 应 用 于 喜马拉雅 地 区 和 喀 喇 昆仑 山西 部 冰 
川 的 研究 。 另 外 ,本 文 用 到 了 第 二 次 冰川 编目 数 
据 、 中 国 行政 区 矢量 数据 。 
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del 遥感 影像 数据 列表 


Tab.1 List ofremote sensing image data 


2.2 研究 方法 

2.2.1 冰川 边界 提取 与 精度 评估 方法 ”很 多 计算 机 
数据 人 处理 软件 提供 了 多 种 自动 或 半自动 提取 冰川 
边界 的 方法 ,但 这 些 方法 的 应 用 在 精度 方面 仍 受 到 
很 大 的 限制 。 比 如 有 阔 值 统计 法 、 非 监督 分 类 和 监 
督 分 类 、 主 成 分 分 析 法 、 面 向 对 象 的 信息 提取 法 、 比 
值 闵 值 法 等 皇宫。 其 中 ,比值 阔 值 法 可 以 在 地 学 处 
理 软 件 ENVI 中 进行 ,原理 简单 .提取 结果 较 准确 ， 
尽管 其 对 冰川 表 研 的 识别 较 差 ,但 仍然 是 目前 应 用 
最 多 且 可 信 度 最 高 的 冰川 提取 分 类 方法 人。 本 研 
究 发 现 ,比值 闵 值 法 受到 影像 质量 、 冰 奔 物 遮挡 、 视 
觉 判 断 误差 等 因素 的 影响 很 大 ,因此 基于 2 次 冰川 
编目 数据 以 及 1989 一 2018 年 遥感 影像 ,借助 Google 
Earth 采用 目 视 解 译 提取 冰川 边界 。 冰 川 边界 精度 
主要 受 传 感 器 和 图 像 配 准 误 差 的 影响 。 本 文 仅 考 
虑 遥感 影像 空间 分 辩 率 造成 的 误差 ” ,可 由 下 式 计 
算得 出 : 


e=NXA (1) 
式 中 : s 为 影像 空间 分 辩 率 造成 的 冰川 面积 误差 
(km?) ;NN 为 冰川 轮廓 的 周 长 (km) ;4 为 半 个 像 元 的 
WK (Landsat TM/ETM+/OLI 2j 15 m)(m)。 结 果 表 
HH, e RAN 39.27 km , 占 2018 年 总 面积 的 2.7% , 
其 他 年 份 的 冰川 边界 精度 也 是 2.7% ,精度 满足 
2.2.2 冰川 流速 误差 评估 方法 ”误差 评估 方法 最 为 
精准 的 是 进行 实验 数据 与 实测 数据 的 对 比 验 证 ,但 
喀 喇 昆仑 出 的 冰川 海拔 高 ` 地 势 险要 , 受 自然 条 件 


序号 日 期 (年 -月 -日 ) 传感器 云 量 序号 日 期 (年 -月 -日 ) 传感器 云 量 
1 1989-11-17 TM 10.18 14 2004-08-14 ETM+ 1.30 
2 1992-11-25 TM 18.74 15 2005-09-02 ETM+ 6.94 
3 1991-08-19 TM 18.87 16 2006-10-07 ETM+ 2:72. 
4 1993-07-07 TM 1.00 17 2007-04-17 ETM+ 2.67 
5 1994-06-24 TM 5:21 18 2008-06-06 ETM+ 1.99 
6 1996-09-01 TM 17.60 19 2009-08-12 ETM+ 1.84 
"i 1997-07-18 TM 0.83 20 2010-10-18 ETM+ 2.77 
8 1998-04-16 TM 0.00 21 2011-08-02 ETM+ 31.64 
9 1999-01-13 TM 0.00 22 2013-10-10 ETM+ 2.17 
10 2000-06-16 TM 2.00 23 2014-06-15 OLI 8.84 
11 2001-05-18 TM 2.00 24 2015-01-09 OLI 11.20 
12 2002-10-12 TM 3.00 25 2016-09-24 OLI 13.38 
13 2003-05-08 ETM+ 2.00 26 2018-04-07 OLI 2.75 


的 限制 ,长 期 以 来 难以 大 面积 的 进行 实地 测量 。 且 
受 影像 配 准 ` 云 雪 质 量 . 冰 面 光滑 、 冰 川 表 磺 物 识 
别 .ETM 影像 修复 、 冰 川 边界 提取 不 确定 性 等 的 影 
响 ,导致 计算 的 不 同时 期 间隔 的 位 移 量 相关 性 不 
好 ,影响 一 定 的 估算 精度 ,可 采用 静止 的 山峰 等 物 
体位 移 为 0 的 非 冰 川 区 为 基准 , 评 佑 估算 精度 *。 
为 此 ,本 文 以 2018 年 为 例 ,根据 遥感 冰川 运动 估算 
的 基本 原理 ,利用 非 冰 川 区 残余 位 移 对 冰川 表面 运 
动 监 测 结果 开展 了 精度 分 析 "" 其 详细 公式 如 下 : 

e, = JMED + SE (2) 
式 中 :eu 为 非 冰 川 区 偏 移 量 的 误差 ;MED 为 非 冰 川 
区 偏 移 量 的 平均 值 ;SE 为 非 冰 川 区 偏 移 量 的 标准 误 
差 , 其 具体 计算 公式 如 下 : 

sE=SIDV (3) 

式 中 :STDV 为 非 冰川 区 偏 移 量 平均 值 的 标准 差 SIN 
为 空间 非 相 关 像 元 个 数 ,其 计算 公式 如 下 : 


(4) 


式 中 :Ni 为 非 冰川 区 像 元 的 总 数 ;PS 为 像 元 分 状 率 
(m);D 为 空间 相关 的 最 大 距离 (m) ,空间 相关 的 最 
大 距离 一 般 是 像 元 分 辨 率 的 20 倍 。 由 此 可 得 非 冰 
川 区 偏 移 量 的 误差 为 1.413 ma', 远 小 于 冰川 的 速 
度 观 测 值 。 

2.2.3 冰川 类 型 识别 ” 喀 喇 昆仑 山北 坡 42 条 冰川 
^H BIARRA i ,使 得 冰川 边界 的 识别 成 为 一 大 
困难 ”。 利 用 Landsat 系 列 影像 和 冰川 编目 数据 绘 
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制 冰 川 末 端 来 决定 冰川 是 否 前 进 LER E SB AR TU 。 
按照 Rankl 等 ”提出 的 定义 :以 2 个 像 元 大 小 为 标 
准 , 判 断 冰 川 末 端 前 进 与 否 。 当 冰川 末端 前 进 超过 
60 m 时 为 前 进 冰川 ;冰川 末端 退缩 超过 60 m 时 为 退 
缩 冰 川 ; 否 则 为 稳定 冰川 外 。 对 于 跃 动 冰 川 的 识 
别 ,可 根据 其 末端 变化 .表面 流速 .表面 特征 (如 裂 
BR ` 冰 面 明显 的 中 矿 裙 皱 、 冰 面 的 破碎 ) 未 端 增 厚 
识别 外。 一 般 而 言 , 跃 动 冰川 冰 流 速度 超出 平静 
期 10 倍 以 上 中 ,往往 造成 末端 突然 前 进 和 冰川 上 部 
的 减 薄 下 。 目 前 ,关于 具体 区 分 跃 动 冰 川 与 非 跃 动 
冰川 还 没有 明确 的 定义 。 并 且 跃 动 与 前 进 冰川 存 
在 一 定 的 区 别 ”。 一 般 认 为 ,未 端 速度 大 于 100.00 
ma 的 前 进 冰 川 .表面 流速 大 于 100.00 m-a” vk 
川 可 归 为 跃 动 。 冰 川 末端 速度 可 根据 2 期 冰川 末端 
长 度 变 化 获取 ,冰川 长 度 的 获取 需要 获取 冰川 中 流 
线 , 冰 川中 流 线 长 度 即 冰川 长 度 。 冰 川中 流 线 可 在 
ArcGIS 中 通过 半自动 提取 获得 。 具 体 方法 是 基于 
冰川 边界 和 DEM 提取 冰川 最 高 点 到 冰川 最 低 点 的 
PDR, 


3 结果 与 分 析 


3.4 冰川 分 类 结果 

基于 1989 一 2018 年 遥感 影像 图 ,根据 冰川 末端 
变化 识别 了 21 条 前 进 冰 川 .20 条 稳定 冰川 1 条 退 
缩 冰 川 。1989 一 2000 年 发 生前 进 的 冰川 有 8 条 ,发 
生 退 缩 的 冰川 有 2 条 ;2000 一 2010 年 发 生前 进 的 冰 
川 有 14 条 ,发 生 退 缩 的 冰川 有 3 条 ;2010 一 2018 年 
发 生前 进 的 冰川 有 13 条 ,退缩 冰川 2 条 ,其 余 冰 川 
保持 稳定 。 在 1989 一 2018 年 期 间 , 存 在 稳定 冰川 16 
条 前进 冰 川 6 条 、 退 缩 冰川 1 条 、 跃 动 冰 川 19 条 ,其 
中 仍 未 停止 前 进 的 冰川 有 4 条 。 末 端 前 进 速 度 大 于 
100.00 ma' 的 冰川 有 11 条 ;前 进 速 度 大 于 200.00 
ma' 的 有 8 条 ;前 进 速度 大 于 300.00 ma' 的 有 8 
条 。 前 进 期 短 则 1a, 长 则 24 a, 发 生前 进 的 时 间 主 
要 集中 在 近 10 a。 跃 动 冰 川 分 布 在 前 进 冰 川 .退缩 
冰川 .稳定 冰川 中 。 根 据 冰川 表面 特征 与 流速 变化 识 
别 了 19 条 跃 动 冰川 。 其 中 ,冰川 末端 发 生前 进 的 路 
动 冰 川 有 12 条 。 图 2 为 不 同类 型 冰川 的 空间 分 布 。 
3.2 冰川 运动 特征 
3.2.1 冰川 表面 运动 速度 室 间 分 布 将 征 ”本 人 研究 基 
于 美国 雪 冰 中 心 的 ITS_LIVE 数 据 对 喀 喇 昆仑 山北 
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图 2 雅 动 冰川 前进 冰川 ,退缩 冰川 和 稳定 冰川 的 空间 分 布 
Fig.2 Spatial distribution of surging glacier, advancing 


glacier, retreating glacier and stable glacier 


坡 42 条 主要 冰川 流速 做 了 全 范围 的 估算 ,得 到 了 
1990 一 2018 年 流速 的 空间 分 布 情况 (图 3)。 以 2018 
年 为 例 ,分析 研究 区 不 同 状 态 的 冰川 运动 特征 。 研 
究 得 到 ; 跃 动 冰川 年 平均 流速 介 于 1.38~38.8 m-a” 
之 间 ,平均 流速 为 14.62 ma'!; 前 进 冰 川 流 速 介 于 
1.11~9.04 m-a :之 间 ,平均 流速 为 4.32 m-a”; FAE VK 
川 流 速 介 于 0.64~7.9 m+ a! Z8] ,平均 流速 为 3.12 
ma'; 退 缩 冰川 流速 为 3.54 ma'。 显 然 , 跃 动 冰 川 
流速 大 于 其 他 3 类 冰川 流速 。 在 空间 上 ,流速 分 布 
符合 冰川 运动 一 般 原理 ,可 在 积累 区 和 冰 舌 上 部 
发 现 明 显 的 高 值 区 域 ,这 些 区 域 多 出 现在 跃 动 冰 川 
主干 或 支部 ,流速 大 小 一 般 在 100.00 m-a ZE, Agr 
苏 盖 提 冰川 南 分 支 ,其 最 大 流速 可 达到 358.33 ma '。 
单条 冰川 平均 流速 分 布 最 大 的 冰川 为 9 号 冰川 
(38.80 m-a”) ,分 布 最 小 的 为 32 号 冰川 (0.64 ma ')。 

冰川 流速 的 空间 分 布 与 高 程 关 系 密 切 。 研 究 
区 冰川 海拔 处 于 3800~8600 m ,我 们 将 海拔 以 每 隔 
400 m 分 段 ,得 到 2018 年 不 同 高 程 区 间 的 冰川 平均 
流速 (图 4)。 结 果 显 示 : 随 着 海拔 升 高 ,冰川 流速 逐 
渐 升 高 后 下 降 。 在 海拔 3800~4200 m 之 间 ,冰川 流 
速 为 3.21 ma', 是 顶部 8200~8600 m 区间 流速 的 2 
倍 ; 接 着 , 随 着 海拔 升 高 ,冰川 流速 从 3.21 ma ' 升 高 
到 54.55 m - a” Jri A Ir va) B 1.58 m - a^ ,在 海拔 
4600-5000 m Z [H] ,流速 达到 最 大 (54.55 mea!) ,是 
冰川 末端 流速 的 10 倍 以 上 ;到 达 5000 m 以 后 ,流速 
逐渐 减 小 ,在 7000 m 以 后 ,冰川 运动 最 慢 。 

冰川 运动 还 与 坡度 、 坡 向 等 有 关 。 坡 度 分 布 在 
0~5° 之 间 的 42 条 冰川 平均 流速 最 大 ,并 且 随 着 坡度 
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图 4 不 同 地 形 参数 的 冰川 流速 
Fig. 4 Glacier velocity with different topographic parameters 
的 增 大 ,流速 逐渐 减 小 。 这 与 Sam EAR 3.2.2 KIADAA ERRE ARS EXC 


9 的 研究 结果 类 似 , 可 能 是 因为 坡度 小 ,利于 冰川 
面积 积累 ,物质 交换 快 ,重力 分 量 大 ,推动 冰 体 快速 
向 下 流动 。 处 于 东 向 的 冰川 流速 最 大 ,例如 处 于 东 
朝向 的 音 苏 盖 提 冰川 ,平均 流速 可 达 21.18 ma ;人 处 
于 西南 向 的 流速 最 小 (17.02 m-a'). 


进行 的 冰川 分 类 ,研究 不 同类 型 的 冰川 流速 年 际 变 
化 特征 。 研 究 区 16 条 稳定 冰川 末端 没有 展现 出 明 
显 的 前 进 和 表面 特征 的 变化 。 在 7 期 4 条 稳定 冰川 
1990—2018 年 中 流 线 流速 剖面 图 (图 5) 中 发 现 , 除 
了 34 号 冰川 有 轻微 波动 外 ,其 他 3 条 冰川 年 际 变化 
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图 5 稳定 冰川 不 同年 际 中 流 线 流速 
Fig. 5 Streamline velocities at different interannual periods of stable glaciers 
均 较 稳定 。 份 流速 保持 一 致 。10 号 冰川 在 距 顶 端 10000~ 


本 文 共识 别 了 6 条 前 进 冰 川 , 分 别 是 4 号 、19 
号 ,20 号 .25 号 ,26 号 ,35 号 冰川 。 绘 制 其 剖面 图 
(图 6) 可 发 现 :冰川 流速 年 际 变化 较 大 ,变化 幅度 在 


20000 m 之 间 年 际 变化 较 大 ,其 中 ,该 条 冰川 在 1995 
All 2015 年 各 发 生 一 次 跃 动 , 跃 动 时 峰值 流速 可 达 
320.00 ma', 而 其 余年 份 流速 在 100.00 m» a` DA 


0~30.00 ma ' 之 间 , 整 体 呈 现 不 稳定 状态 。 如 图 7 
展示 了 19 号 和 26 号 冰川 末端 前 进 状 况 。19 号 冰川 
在 2013 一 2014 年 前 进 距离 最 大 , 共 94.82 m。26 号 
冰川 在 1994 一 2007 年 发 生前 进 , 前 进 速 度 可 达 
24.46 ma'!, 其 中 ,1999 一 2000 年 前 进 距离 最 大 ,前 
进 速度 可 达 71.80 mea’. 

退缩 冰川 (图 8)12 号 冰川 在 1993—1996 年 之 
间 发 生 退 缩 ,可 以 发 现 其 流速 波动 变化 ,变化 幅度 
在 0~12.00 ma 之 间 , 整 体 呈 不 稳定 。 

跃 动 冰川 有 2 种 行为 特征 。 一 种 是 末端 没有 发 
生前 进 的 ,一 种 是 末端 发 生前 进 的 中。 如 图 9~10 所 
示 : 部 分 末端 稳定 的 冰川 在 积累 区 或 消融 区 表现 出 
明显 的 跃 动 。 例 如 0 号 和 10 号 冰川 。0 号 冰川 在 距 
冰川 顶端 1000~2000 m 之 间 的 消融 区 位 置 流 速 大 于 
100.00 ma ,并 在 该 区 出 现年 际 差异 ,其 表现 为 
1990 年 和 2010 年 流速 大 于 其 他 年 份 流速 ,而 其 他 年 


下 。2 号 冰川 和 17 号 冰川 在 积累 区 年 际 变化 较 大 。 
2 号 冰川 距 顶端 5000~15000 mm 之 间 年 际 变化 较 大 ， 
2000 年 流速 达到 了 最 低 ,1990.2015 年 和 2018 年 流 
速 在 距 顶 端 7500~12000 m 之 间 都 超过 100.00 m-a”, 
说 明 该 条 冰川 至 少 发 生 了 3 次 跃 动 。 不 同 冰川 不 同 
位 置 年 际 流速 变化 均 有 差异 。17 号 冰川 年 际 变化 
较 复杂 ,该 条 冰川 在 距 顶 端 5000~25000 m 之 间 年 际 
变化 较 大 ,而 在 距 顶 端 253000 m 到 冰川 末端 年 际 变 
化 稳定 , 且 不 同年 份 的 流速 时 空 变化 一 致 。 其 中 ， 
1995 ^E fll 2015 4F vk) || TERR Jit ig 5000-15000 m zz fH] 
发 生 了 2 次 跃 动 ;1990 年 和 2010 年 冰川 在 距 顶端 
15000-25000 m 之 间 也 各 发 生 一 次 跃 动 。5 号 和 9 
号 冰川 年 际 变化 较 小 。27 号 .39 号 .40 号 冰川 年 际 
变化 较 大 ,但 流速 峰值 出 现 的 位 置 与 往年 大 致 相 
同 。39 号 冰川 在 2015 年 和 2018 年 发 生 过 2 次 跃 动 ， 
跃 动 周期 小 , 跃 动 前 端 基本 未 变 ,但 冰川 表面 特征 
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图 6 前 进 冰 川 不 同年 际 中 流 线 流速 


Fig.6 Streamline velocities at different interannual periods of advancing glaciers 
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7 典型 冰川 前 进 现象 


Fig.7 Typical glacial advance phenomena 


发 生变 化 ,运动 将 冰 与 表 矿 的 混合 物 向 下 输送 沾 。 ”类 冰川 中 流 线 流速 的 年 际 变 化 不 大 ,在 0~20.00 ma 
末端 发 生前 进 的 冰川 ,大 多 是 中 规模 的 冰川 ,例如 7 ”之 间 变 化 ,整体 呈现 不 稳定 现象 。 尽 管 这 类 冰川 表 
号 .18 号 . 22 号 . 28 号 ,29 号 30 号 .33 号 .37 号 。 这  ” 面 流 速 较 小 ,但 其 末端 运动 速度 较 大 ,会 发 生 几 百 
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图 8 退缩 冰川 12 号 冰川 不 同年 际 中 流 线 流速 
Fig. 8 Streamline velocities at different interannual periods 


of retreating glacier No.12 


米 的 前 进 。 如 图 11 展 示 了 典型 跃 动 冰 川 在 1989 一 
2018 年 的 末端 加 速 情况 。 如 33 号 冰川 , 跃 动 启动 时 
间 在 1993 年 ,结束 在 2018 年 , 跃 动 期 可 达 25 a, 未 端 
在 2018 年 至 往 后 仍 有 继续 前 进 的 趋势 。37 号 冰川 
跃 动 启动 时 间 在 1993 年 ,结束 在 2010 年 , 跃 动 期 可 
达 27 a, 峰 值 速度 出 现在 2009 一 2010 年 之 间 , 流 速 
达到 586.84 ma'。 


4 讨论 


4.1 冰川 规模 与 流速 

从 表 2 可 得 :研究 区 面积 在 1.71~359.05 km? 之 
间 变 化 。 跃 动 冰 川 主要 发 生 在 面积 大 于 5 km 以 上 
的 冰川 ,而 稳定 冰川 多 集中 在 面积 在 2.0~5.0 km 之 
间 的 冰川 。 面 积 大 于 100 km 的 冰川 有 4 条 ,平均 流 
速 分 别 为 32.82 m*a/,38.80 ma 21.18 m-a, 
16.95 mea’; 面积 分 布 在 50~100 km 之 间 的 冰川 有 3 
条 ,平均 流速 分 别 为 25.30 m+ a! 13.73 ma'、35.05 
ma!。 这 7 条 冰川 的 流速 都 大 于 面积 <50 km 的 冰 
川 的 流速 。 跃 动 冰 川 的 长 度 一 般 要 大 于 其 他 3 类 冰 
川 , 即 冰川 跃 动 更 容易 发 生 在 长 度 更 长 的 冰川 上 ， 
而 跃 动 冰川 的 流速 要 大 于 稳定 冰川 流速 。 冰 川 规 
模 与 表面 流速 具有 一 定 的 关系 。 面 积 大 的 冰川 有 
大 面积 的 积累 区 ,积累 区 有 丰富 的 冰川 物质 ,促使 
冰川 向 前 运动 。 因 此 ,面积 大 的 冰川 流速 较 大 。 
4.2 气候 对 冰川 流速 的 控制 

井 哲 帆 等 “通过 对 冰川 形态 物质 平衡 \ 气 候 条 
TE .冰川 面积 和 厚度 等 影响 冰川 运动 的 因子 进行 通 
径 分析 ,得 到 影响 冰川 运动 速度 的 因子 相互 影响 ， 
其 中 直接 影响 冰川 运动 速度 的 因子 是 冰川 形态 ,而 


气候 因素 则 是 非常 重要 的 间接 因素 ,通过 气候 的 变 
化 影响 物质 平衡 ,从 而 作用 于 冰川 形态 以 及 一 系列 
冰 内 结构 与 应 力 影 响 冰 川 运动 速度 。 不 同 的 地 理 
位 置 ,不 同 的 冰川 ,其 影响 运动 速度 的 因素 是 不 相 
同 的 。 气 候 因 素 是 控制 因子 ,气候 的 变化 控制 着 冰 
川 的 运动 状态 。 

喀 喇 昆仑 山北 坡 所 处 的 地 理 位 置 特殊 ,受到 西 
风 环 流 和 印度 洋 西 南 季 风 影 响 , 冬 、 夏 两 季 会 有 大 
量 的 降水 ,降水 对 冰川 物质 的 积累 有 重要 的 有 影响。 
我 们 通过 许 艾 文 等 所 的 研究 得 到 ,1960 一 2013 年 各 
站 点 降水 量 逐 渐 增 加 ,年 均 气 温 与 夏季 气温 都 呈现 
显著 上 升 趋势 。 虽 然 整 体 趋势 相同 ,但 变化 幅度 也 
有 差异 。 冰 川 发 育 的 条 件 相 比 来 说 会 影响 冰川 的 
发 育 及 稳定 状态 ,如 果 破 坏 了 这 种 稳定 平衡 的 发 育 
ARTE ,就 极 有 可 能 造成 冰川 的 不 稳定 。 因 此 , 喀 喇 
昆仑 山 的 这 种 气候 变化 差异 是 造成 一 部 分 地 区 冰 
川 不 稳定 .不 平衡 ,从 而 发 生前 进 或 跃 动 的 主导 因 
素 。 冰 川 运动 受气 候 控 制 , 最 直观 的 现象 就 是 冰川 
流速 的 季节 性 变化 。 当 降水 增加 ,冰川 积累 增加 ， 
厚度 增加 ,推进 冰川 运动 。 当 气温 升 高 时 ,冰川 融 
水 增加 ,促进 底部 滑动 ,影响 冰川 流速 。 研 究 区 冰 
川 流 速 在 1989 一 2018 年 呈 增 大 趋 热 。 这 与 研究 区 
增加 的 降水 与 气温 关系 密 不 可 分 。 
4.3 冰川 跃 动机 制 

冰川 跃 动 是 一 种 复杂 的 活动 ,至今 未 有 对 其 机 
制 的 统一 解释 “” 。 现 有 研究 对 喀 喇 昆仑 山 冰 川 跃 
动 的 原因 有 不 同 看 法 ,认为 喀 喇 昆仑 山 冰 川 跃 动 受 
热力 学 和 水 文 条 件 控制 。 井 哲 帆 等 “指出 气候 因素 
是 控制 因子 ,气候 的 变化 控制 着 冰川 的 运动 状态 。 
但 从 Rankl 等 中 和 Quincey 等 的 研究 结果 可 以 得 
到 , 喀 喇 昆仑 山 存 在 较 多 的 跃 动 型 冰川 ,这 种 冰川 
本 身 也 会 引起 局 部 冰川 条 件 的 改变 。 冰 川 跃 动 也 
有 可 能 是 冰川 内 部 或 底部 岩石 受到 冰 体 挤 压 ,以 及 
水 文 条 件 引 起 冰 温 变化 的 综合 作用 。 当 冰 面 坡度 
较 大 时 ,发 生 朋 裂 , 向 下 滚 落 ,造成 整个 冰 体 不 稳 
定 ,从 而 快速 下 请 造成 冰 流 速度 增 大 而 发 生 跃 动 。 


5 结论 


本 文 利 用 ITS_LIVE .Landsat 数 据 提 取 了 研究 区 
42 条 冰川 的 表面 流速 ,并 将 末端 运动 与 表面 运动 特 
征 结 合 起 来 ,对 比分 析 了 不 同 运 动 类 型 的 冰川 表面 
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图 9 跃 动 冰川 不 同年 际 中 的 流 线 流速 (冰川 序号 :0、2、5、7、9、10、17、18、21、22、27、28) 


Fig.9 Streamline velocities at different interannual periods of surging glaciers 
(glacier numbers: 0.2.5.7.9.10.17.18.21.22.27.28) 
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图 10 跃 动 冰 川 不 同年 际 中 的 流 线 流速 (冰川 序号 :29、30、31、37、39、40) 


Fig. 10 Streamline velocities at different interannual periods of surging glaciers (glacier numbers: 29、30、31、37、39、40) 
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图 11 典型 冰川 跃 动 现象 


Fig. 11 Typical glacier surging phenomenon 
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表 2 选取 的 典型 冰川 基本 信息 表 


Tab.2 Basic information for the selected typical glaciers 


冰川 编号 冰川 编码 坡 向 面积 km” ”长 度 /km 平均 高 程 /m 平均 流速 /ma 冰川 类 型 
0 5Y654C0163 NE 107.74 21.98 5670.75 32.82 跃 动 
1 5Y654C0167 E 32.25 13.83 5587.68 6.38 稳定 
2 5Y654C0157 NE 87.83 26.91 5358.23 25.30 ERZ) 
3 5Y654C0136 N 1.71 3.22 5488.77 1.66 稳定 
4 5Y654C0045 NW 6.14 5.89 5793.50 3.53 前 进 
5 5Y654C0145 N 83.82 28.83 5519.46 13.73 跃 动 
6 5Y654C0135 NE 2.03 3.54 5572.33 2.74 稳定 
7 5Y654C0132 NE 10.00 8.42 5574.16 3.32 ERZ) 
8 5Y654C0130 NE 6.16 8.72 5537.67 2.70 稳定 
9 5Y654C0128 N 110.35 26.42 5617.93 38.80 ERZ) 
10 5Y654C0116 NE 94.43 20.12 5664.67 35.05 跃 动 
11 5Y654C0092 N 39.53 14.53 5699.06 7.86 TAXE 
12 5Y654C0087 NE 6.93 6.81 5807.98 3.54 退缩 
13 5Y654D0123 NE 2.34 4.29 5534.80 1.78 稳定 
14 5Y654D0122 N 2.72 3.19 5458.10 2.45 稳定 
15 5Y654D0119 NE 7.43 5.44 5507.69 2.64 稳定 
16 5Y654D0112 NE 9.64 8.15 5443.55 7.30 TAXE 
17 5Y654D0053 E 359.05 43.16 5377.05 21.18 ERZ 
18 5Y654D0097 E 19.80 11.37 5660.93 8.24 跃 动 
19 5Y654D0105 NE 10.16 6.58 5529.76 5.05 前 进 

20 5Y654D0102 NE 5.92 5.98 5597.82 2.57 前 进 
21 5Y654D0096 N 27.55 15.10 5418.99 9.69 跃 动 
22 5Y654D0064 8.17 6.85 5756.31 1.38 跃 动 
23 5Y654D0092 NE 241 4.86 5686.01 2.57 稳定 
24 5Y654D0089 NE 3.78 5.05 5596.15 2.32 稳定 
25 5Y654D0066 S 1.76 3.76 5712.78 1.11 前 进 
26 5Y654D0087 NE 12.20 6.46 5519.86 9.04 前 进 
27 5Y654D0048 NE 182.54 27.91 5450.89 16.95 ERZ) 
28 5Y654D0048 E 8.54 8.08 5509.20 2.56 跃 动 
29 5Y654D0077 NE 11.92 7.65 5524.96 2.88 跃 动 
30 5Y654D0068 SE 10.09 10.64 5792.27 4.39 跃 动 
31 5Y654D0073 N 3.26 5.15 5569.96 1.29 稳定 
32 5Y654D0072 NE 241 4.57 5503.74 0.64 稳定 
33 5Y654D0033 NW 11.14 7.03 5826.72 5.42 跃 动 
34 5Y654D0032 NW 4.50 5.54 5945.17 1.25 稳定 
35 5Y654D0027 N 3.37 5.53 5660.38 1.67 前 进 
36 5Y654D0026 NW 2.13 4.46 5849.63 0.95 稳定 
37 5Y654D0023 NW 3.43 6.34 5475.77 1.92 跃 动 
38 5Y654D0022 N 1.73 3.87 5438.58 1.45 稳定 
39 5Y654D0016 NE 74.68 21.31 5778.15 31.95 跃 动 
40 5Y654D0016 16.44 12.74 5721.73 7.51 跃 动 
41 无 N 38.03 15.37 5511.19 11.29 稳定 
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流速 时 空 变化 。 得 到 如 下 结论 : 

(1) 1989 一 2018 年 ,人 研究 区 42 条 冰川 中 ,存在 
稳定 冰川 16 条 、 前 进 冰川 6 条 、 退 缩 冰川 1 条 、 跃 动 
冰川 19 条 。 在 此 期 间 , 仍 未 停止 前 进 的 冰川 有 4 
条 。 末 端 前 进 速度 大 于 100.00 ma' 的 冰川 有 11 
条 ;前 进 速度 大 于 200.00 ma'! 的 有 8 条 ;前 进 速 度 
大 于 300.00 ma' 的 有 8 条 。 前 进 期 短 则 1 a, 长 则 
24 a, 发 生前 进 的 时 间 主 要 集中 在 近 10 ao 

(2) 空间 上 ,流速 分 布 符合 冰川 运动 一 般 原 理 ， 
可 在 积累 区 和 冰 舌 上 部 发 现 明显 的 高 值 区 域 ,这 些 
区 域 多 出 现在 跃 动 冰川 主干 或 支部 ,流速 大 小 一 般 
在 100.00 ma 之 上 ,如 在 音 苏 盖 提 冰川 南 分 支 ,其 
最 大 流速 可 达到 358.33 ma'。 单 条 冰川 平均 流速 
分 布 最 大 的 冰川 为 9 号 冰川 (38.80 m-a”) ,分 布 最 小 
的 为 32 号 冰川 (0.64 ma). 

(3) 整体 上 ,冰川 与 地 形 要 素 的 关系 密切 。 流 
速 在 海拔 4600~5000 m 之 间 达 到 最 大 (54.55 m-a"), 
是 冰川 末端 流速 的 10 信 以上; 分布 在 坡度 0~5° 之 间 
的 42 条 冰川 平均 流速 最 大 ,并 且 随 着 坡度 的 增 大 ， 
流速 逐渐 减 小 ;处 于 东 向 的 冰川 流速 最 大 ,例如 处 
于 东 朝 向 的 音 苏 盖 提 冰 川 ,平均 流速 可 达 21.18 me 
al ,处 于 西南 向 的 流速 最 小 (17.02 ma). 

(4) 不 同类 型 的 冰川 中 流 线 流速 年 际 变化 不 
同 。 稳 定 冰川 流速 年 际 变化 较 稳 定 ,不 同年 份 相同 
位 置 的 流速 值 较 一 致 ; 而 对 于 前 进 冰川 和 退缩 冰川 
而 言 , 年 际 流速 波动 均 较 大 ; 跃 动 冰川 不 同位 置 年 
际 流速 变化 不 同 。 末 端 发 生前 进 的 跃 动 冰川 表面 流 
速 变 化 不 大 ,但 末端 变化 速度 却 可 超过 100.00 m-a”, 
甚至 可 达到 500.00 m-a", 

(5) 长 度 更 长 .面积 更 大 的 冰川 易 发 生 跃 动 。 
气候 、 雪 月 和 地 形 也 是 影响 冰川 运动 速度 变 大 的 主 
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Abstract: The surface velocity of 42 glaciers on the northern slope of K2 (Mount Qogir), Karakoram was 
extracted from ITS LIVE and Landsat data. Combined with data regarding the characteristics of end and surface 
movements, the spatial and temporal variations of the surface velocity of glaciers with different movement types 
were analyzed and compared. The results show that (1) from 1989 to 2018, among the 42 glaciers in the study 
area, we observed 16 stable, 6 advancing, 1 retreating, and 19 surging glaciers. The velocity distribution in space 
is consistent with the general principle of glacier movement, and obvious high-value areas can be found in the 
accumulation area and upper region of the glacier tongue; these areas are mostly in the trunk or branch of the 
active glacier, and the velocity is generally 2100.00 m-a '. For example, in the southern branch of the Yengisogat 
glacier, the maximum velocity can reach 358.33 m-a '. (2) The glacier velocity distribution is closely related to 
topographical factors. The maximum flow velocity (54.55 m +a ') is between 4600 m and 5000 m, which is >10 
times that of the end of the glacier. The average velocity of the 42 glaciers in the slope angle range of 0? to 5? is 
the highest, and the velocity decreases with the increase of the slope angle. Glaciers flowing toward the east have 
the greatest flow rate, and those that flow toward the southwest have the smallest. (3) The interannual variation of 
the velocity of the stable glacier is small, and the velocity values at the same location are consistent across 
different years. For advancing and retreating glaciers, the fluctuation in interannual velocity is large. The annual 
velocity variation of the surging glacier is found to be location dependent. The surface velocity of the glacier with 
a moving tip shows little change, but the change velocity of the glacier tip can be 7100.00 m-a ' and even up to 
500.00 m :a '. (4) Glaciers with longer lengths and larger areas are prone to surge. Climate, avalanches, and 
topography also play a major role in determining glacier speed. 

Key words: glacier; the northern slopes of the Karakoram Mountains; ITS LIVE data; Landsat data; surface 
velocity; the spatial and temporal variations 


